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Summary 

During research period the following models of the VASIMR helicon discharge have 
been further developed and applied to analyze the on-going VX-10 ASPL experiment: 

A) ID semi-analytical model for a mixed-collisional propellant flow 

B) OD power and balance model for the whole helicon discharge 

In this particular research period we have concentrated on the MW-level performance of 
the VASIMR helicon source. Favorable high-power scaling and reduced ionization costs 
were obtained, and presented at the VASIMR NASA review in the Fall ’02. This Grant is 
continuation of the previous NAG9-1224 award. The research results are summarized in 
14 publications; they were presented as 20+ talks at the major International Conferences 
and scientific seminars at the leading Academic and Research Institutions. The reported 
results allowed helicon discharge characterization, understanding of the several 
experimental observations, and helped to make predictions and propose structural 
modifications for the advanced VASIMR helicon source operation. 
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Motivation 


The integral performance of the high, variable specific impulse I sp thruster being 
developed under the VASIMIR project [1] by ASPL/NASA depends on the efficiency of 
the helicon plasma source. Current ASPL experimental data shows stable helicon source 
[1-2] performance for various gases. However, there is new evidence that current models 
and theories of the helicon discharge are incomplete. For instance, RF-power absorption 
may have several resonances; performance with mixtures of gases is different from pure 
hydrogen or helium; retarded probe analyzers show a large population of energetic ions 
with non-Maxwelian distribution; neutral pressure jumps by a factor of 4 during 
discharge; heat deposition onto inner wall has a non-monotonic spatial profile. There is 
an indication that a fraction of RF-power goes into ions, and that charge-exchange 
collision process and residual vacuum tank pressure influence the dynamics of the helicon 
discharge. Because experimental diagnostics are limited new theoretical-numerical 
models are required to understand VX experiment helicon operation. 

The need for the development of various models comes from the analysis of the typical 
VASIMR discharge conditions, listed below 


Parameters of the VASIMR plasma 

Dimensions of the VX-10 helicon source are about 10x100 cm [2]; 

Magnetic field is non-uniform, varies in the 10-600 Gs range; 

Neutral gas pressure varies in the 1-600 mTorr range; 

Neutral gas density alters from 5xl0 15 cm' 3 to lxlO 13 cm' 3 ; 

Plasma density is about l-2xl0 12 cm' 3 
Electron temperature is about 5-10 eV; 

Ion exhaust velocity is about 40 km/sec, corresponding to flow energy ~16eV; 
Plasma flow is sub-sonic or super-sonic; 

Heating power is in 3 KW-10KW range (and up to MW-level); 

RF frequency is 13.6 MHz (and up); 

Mass flow rate is 10-1000 seem. 
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Introduction 


VASIMR is a unique electrode-free electric propulsion system [1-4], which has a potential of 
achieving specific impulse on the order of 30Ksec with 50% energy efficiency. Fuel efficiency 
depends entirely on that of a plasma source. Presently the VASIMR project considers the helicon 
discharge plasma source [5,8] a primary candidate. 


A typical VASIMR thruster combines three distinct areas i) plasma source, ii) plasma RF-heating 
cell and iii) magnetic nozzle, as shown in the Fig.l. The most probable propellant is H 2 to 
achieve maximum specific impulse Is P at a given mass fueling rate. There are several important 
collisional processes taking place in such a thruster. Firstly, there are a few inelastic channels to 
break incoming molecular gas into neutral atoms, ionization and excitations of neutrals, and wall 
interactions [6-7]. 



Figure 1 - Plasma facing part of a typical VASIMR engine 


Next, there are a number of elastic collisions of potential importance: Coulomb interaction, 
resonance charge-exchange, etc [9-10]. The typical Knudsen number for the different elementary 
processes varies in a 0.1-10 range, which indicates that the short mean free path limit [6] fails, 
and a pure kinetic analysis is required [11-12]. There are a few approaches to that. 

Comparison of particle-in-cell and finite difference methods shows that the Fokker-Planck 
approach is computationally more efficient [22]. Particle methods suffer from the intrinsic 
statistical noise, which occurs because of the limited number of model particles per cell, and it is 
always 6-8 orders of magnitude larger than the “natural physical” noise. Finite-difference 
methods are free of the statistical fluctuations in the first place, and allow studying both coarse 
and fine effects. From the mathematical point of view numerical particle solution is always non- 
converged (oscillations recede as an inverse square root of the number of trajectories), while grid 
methods may give a converged solution up to the computer accuracy ~10' 14 . 

Next, finite numbers of particles per cell in the numerical simulation make it impossible to 
resolve accurately enough an important energetic tail of the electron and ion distribution 
functions, because the tail is normally exponential function of energy. Practical limit of the 
resolution (for 10-100K particles per cell) is 3-4 thermal velocities. Time averaging helps for 
steady-state regimes only. Thus, non-linear time-dependent study of the transient regimes is 
impossible in principle. 

Grid method with a non-uniform mesh in the velocity space allows equally accurate resolution of 
the cold core and very energetic tail (100 thermal velocities) of the distribution function. These 
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tails are important for plasma diagnostics with probes and neutral beams; they also determine 
plasma thermal conductivity. 

Lastly, in the regimes of interest, regions of sharp spatial gradients of plasma parameters are 
expected. The universal tool to study such fronts is provided by adaptive grids in space [23]. 
While the pure kinetic approach is an ultimate goal, we are also developing simplified 0-D 
plasma chemistry model, which allows us to obtain crucial parameters of the helicon discharge: 
plasma and neutral species composition, temperatures of all plasma components, energy/particle 
fluxes out of the plasma source and heat fluxes onto wall. This model is a continuous 
development of our original model [16]. The latest additions include possibility of Lyman-a 
radiation reflection and recapture, as well as it artificial reduction to study the energy efficiency 
of the VASIMR thruster. 

This simple 0-D model uses critical information about propellant flow in the system of pipes as 
given by our 1-D semi-analytical hybrid model, which was also developed under this award. 
New fully kinetic 2V2D BGK-type model for a rarified gas flow independently verifies its results 
in a system of connected channels. Information about plasma flow is taken from our fully kinetic 
Fokker-Planck model that was described in the early publications [16-17]. These results are 
compared to recent RPA data [21] for the ion distribution function obtained for the VX-10 
experimental setup. 

All these numerical models have been applied to: 

a) understand operation of the helicon plasma source in the VASIMR VX-3 and VX-10 
experiments; 

b) make predictions for the future space VF-10/25 thruster [4,5], and 

c) propose practically feasible suggestions on how to boost fuel and energy efficiencies of the 
critical elements of the plasma thrusters, based on the VASIMR technology. 

Current status of VX-10 [5] experiment shows about 40% fuel efficiency for hydrogen, 
deuterium and helium discharges when they operate in a high-density mode. It’s not clear what 
determines transition to this important mode. We use a few numerical models to study 
operational regimes of the helicon plasma source. 

Energy and particle balances are studied using 0-dimentional model of gas discharge [16]. The 
following processes are accounted for hydrogenic discharge: molecule and molecular ion 
dissociation, neutral ion ionization and excitation, charge-exchange, wall conversion, formation 
of heavy ions, etc. Model predicts plasma and gas species densities and temperatures to be very 
close to that measured in the experiment. The model is used to assess effects of gas baffles, 
vacuum tank, gas pre-heating, altering dimensions of the quartz tube, etc. on the hydrogen and 
helium discharges performance. 

Finally we discuss fully kinetic 1D2V and Fokker-Planck-Boltzmann models of the gas 
[11,18,24] and plasma flows [12, 25-26]. They account for spatial variation of ambipolar 
potential and external magnetic field and basic elastic and inelastic collisions. Kinetic 
simulations show that the distribution functions of both gas and plasma species differ from the 
Maxwellian due to external heating, non-local transport and trapping in the magnetic field. IDF is 
in concert with recent RPA measurements [21]. Our numerical results indicate that specific 
impulse on the order of 3-4Ksec is achievable by the helicon source along, and 8-10Ksec with 
additional ICRF-heating of the ions for the 3-10kW input electrical power. 
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A) ID semi-empirical model for propellant flow in pipes 


Our initial numerical and theoretical studies [16-17] have indicated that the degree of gas 
ionization in the helicon source is low, on the order of 1%. This fact allows viewing propellant 
flow in the gas feeding - plasma source system ([19], shown in the Fig.2) as entirely independent 
of that of plasmas. 
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Figure 2 - Gas feed - quartz tube of VX-10 VASIMR configuration 


Simple estimate show that gas flow in the system will be strictly laminar. Indeed, from the 
experimental measurements: pressure drops along the system from PI=100mtorr in the gas inlet 
to PVC~lmtorr in the vacuum chamber, and mass flow rate |j.=80sccm we find that the 
characteristic Reynolds number of the flow 


Re= V d / v = 20 « Rct =2300 (1) 

is much smaller than the critical value [14], Re- , which marks the onset of a turbulent regime. 

On another hand simple consideration shows that the Knudsen number [13] of the propellant 
flow 


Kn=2A/d« 0.5-2 (2) 

Hence, we can view the gas flow to be a composition of the viscous and free molecular laminar 
flows (as shown in the Fig.3 below). 



Figure 3 - Scheme of hybrid model 
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Hybrid Poisseuille - Knudsen model for mixed-collisional gas flow can be easily derived on the 
basis of two classical results. For the viscous flow one can apply Poisseuille result [14] for a 
pipe flow: mass throughput in a pipe is proportional to the pressure gradient and to 4-th power of 
the pipe’s radius R: 

[Ip = A(T) R 4 (z) dP/dz (3) 

where parameter A(T) contains viscosity, etc. 

For a similar free molecular flow one can use analogous result, obtained by Knudsen model [13]: 

p K = B(T)R 3 (z) dP/dz (4) 

Now, taking into account that mass throughput is constant along axial direction, and that 
asymptotic expressions (3)-(4) should be valid for the limits Kn-»0 and Kn— *», respectively, we 
may find the following hybrid expression: 

(J. =/jtr 2 p V z dr = const =p P + ((J. K - Up ) [ 1- exp {-Kn} ] = C(T,P,z) dP/dz (5) 

where C(T,P,z) is a new composite function. By inverting Eq. (5) we obtain the following finite- 
difference expression: 


AP= |x / c(T,P,z) Az (6), 

which can be integrated numerically, say from the gas injector (P=PI) up to the end of the quartz 
tube of the helicon source. Results for He are shown in Fig.4. 



Figure 4 - Profiles of helium gas flow parameters for the VX-10 geometry 
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As one can see from the bottom chart in the Fig.4 agreement with experimental data for the 
pressure drop along the system are within 1% of accuracy for constant helium flow rate 
|i=200sccm. Note that axial temperature profile is a fit to the experimental measurements of the 
quartz tube temperature. We assume that as neutral gas has the same temperature profile due to 
energy exchange with the wall. 

Interestingly enough 80% of the pressure drop occur in the short (8cm) and narrow (1.5cm in 
diameter), 20% in the quartz tube (Fig.2), and under 1% in the rest sections of the system 
(magnet bores, etc.). One can see that our original assumption about mixed viscous-free 
molecular flow is verified a posteriori (as de-facto Knudsen number, Kn, is in the 0.5-5 range). 
The helium gas flow experiences transition viscous — ^molecular — > viscous — > molecular 3, and 
accelerates 2 times. What is the most important for our immediate goal is that mean gas flow is 
very subsonic with average Mach number <M> *0.03-0.05. 

To crosscheck we benchmark hybrid model against experimental data for molecular hydrogen 
gas flow at p=100sccm in the same VX-3 [16] configuration. Corresponding results are close to 
those in Fig.4. 

The major results remain the same: gas flow is very viscous and subsonic, gas density in the 
center of the quartz tube (where helical antenna is located) is on the order of N«3-5xl0 14 cm‘ 3 . 


Comparison of old and new quartz tubes 

Recent results from the VX-10 equipped with a longer quartz tube show two-fold increase of the 
plasma density in the discharge [19]. We apply ID model to study gas flow in the original and 
new tubes. Their geometrical parameters are as follows. Old tube is 5cm in diameter and 70cm 


long. New tube is lm long and has slightly larger 
diameter of 5.5cm. 

Though, dimensions of the two tubes are 
different, surprisingly enough they have the same 
resistance to the hydrogen gas flow. As one can 
see from the Fig.5 both pressure and density 
drops along the tubes coincide with in a few 
percent. Also in the experiment the RF-antenna 
position was fixed at the same position closer to 
the left side of the tube, and the input power 
stayed the same at about 3 kW. 

Thus, the only difference of the two discharges is 
in the gas density at the location of RF-energy 
deposition source. With the longer tube the 
neutral gas density gradient is shallower, energy 
is released to a more dense gas (and plasma too). 
Because gas density is higher, we can guess that 
the ionization may be more efficient, and the 
discharge will yield more plasma for the 
equivalent amount of energy transmitted into the 
electron species. 




Figure 5 - gas profiles in old and new quartz 
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Use of conical tubes instead of cylindrical 

If these arguments are correct, we can try finding simple ways to keep neutral gas density higher 
in the quartz tube volume. One of the possibilities to keep gas density flatter as it accelerates 
axially - is to use a quartz tube of varying cross-section. Because the helicon RF-antenna has 
fixed dimensions we propose using conical shape quartz tubes with average radii approximately 
equal to that of the internal RF-antenna radius. In this case two-way use of the new tubes 
becomes possible automatically. 


A n, 1 0 14 cm' 3 



Figure 6 - gas density in conic tubes 


In Fig.6 we present calculated gas density profiles for a set of three tubes with the same length of 
lm and the same mean internal radius of 5.5cm. Their small-large radii are as follows: 6x5cm, 
6.5x4.5cm and 7x4cm. For comparison we present data for a straight pipe (single linear curve in 
Fig.6). Each pair of curves - one is folded in and another out - correspond to the same conic- 
shape tube oriented in two different ways (gas expansion or compression). As one can see simple 
change of the tube axial orientation allows nearly 100% variation of the gas density in the middle 
of the pipe. On another hand, when the tube’s cross-section is reducing along the flow direction, 
gas density has much more flat profile. In this case we may expect more efficient gas ionization 
by electron impact, and higher overall fuel efficiency of the VASIMR thruster. Note that if we’ll 
take into account gas loss due to ionization, the reduction of tube radius has to be large to 
compensate it. Shaping of the external magnetic field will be obviously required (magnetic 
choke) to avoid hot magnetized plasma striking the internal surface of the tube, causing energy 
loses and unwanted damage. 
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B) Mass and energy balance model for VASIMR plasma source 


The helicon has proven to be a robust and efficient plasma source [8]. The helicon source 
consists of a dielectric (quartz) tube embraced by the helicon antenna, which launches 
electromagnetic waves in the plasma. Electrons are heated through possibly collisional or Landau 
damping of the waves in the plasma. This process is not included self-consistently here; it is a 
subject of a separate study. In this model we assume that a certain power, W, is transmitted to the 
electrons. The electron temperature and density in the source are 6eV and 2xl0 12 cm' 3 , 
respectively, as measured in the ASPL experiments [5, 19]. 


RF electron heating 



Figure 16 - Model of a helicon discharge 


Let’s accept the following physical model of a helicon plasma source, as shown in Fig. 16. Next, 
we consider the accepted models for the main physical processes taking place in the discharge. 

Electron and ion transport . The electrons are magnetized by a relatively strong magnetic field of 
0.1-0.6T. The electron gyroradius is about 10-2cm, ion gyroradius - 0.2cm, which is much 
smaller than 10cm - tube’s diameter. The gyroradius becomes huge in the nozzle. Cross-field 
transport is weak because of the low plasma density. The plasma flux onto outer wall is 
negligible. Electrons strictly follow field lines, which are parallel to the cylinder axis. At one of 
the tube’s ends there is a floating potential wall. The floating potential value for hydrogen 
(deuterium) plasma is 2.8(3) Te. Thus, the majority of the electrons are reflected back. There is a 
small residual particle/energy flux onto the end plate. We anticipate that the continuous puffing 
of H 2 through the inlet will keep the plasma density low near the end plate. Therefore, we neglect 
the heat and mass flow onto the floating end plate. At the open end of the quartz tube the plasma 
is leaving the system at a fraction of sound speed. 

Plasma exhaust velocity . The Debye length of the plasma in typical Helicon hydrogen discharge, 
X D , is on the order of 10' 3 cm, while the tube dimension is on the order of lm. Therefore, plasma 
in the volume is quasineutral and the plasma flow is ambipolar. We assume that the ion’s exhaust 
velocity is sub-sonic with a fraction of the ion sound speed, C S i=(2T e /M i ) 0 ' 5 . The electron density 
is equal to the combined density of all ion species (H\ H 2 + , H 3 + here). The electron exhaust 
velocity is such to automatically maintain the ambipolarity of plasma V e =(n H+ CsH+ + n H2+ C S H2+ + 
nH3+CsH3+V (hh+ + n H2+ + n H 3+). Neutral exhaust velocity is assumed to be sub-sonic as well with 
C S n=0.03 (2Tn/M n ) 0 ' 5 . Because the neutral temperature is significantly lower than the electron 
temperature, ions exhaust the helicon source at a much higher pace than the cold neutrals. 


9 




Energy exhaust and wall fluxes . Sonic neutral exhaust carries along energy flux: 


F F 1 2 

Q £ = oT (CkT + -MC 5 ) = cr(C + 1) nkT V T 


( 7 ) 


where a is a quartz tube cross-section, VT is thermal velocity, C is specific heat. Wall neutral 
flux requires separate attention. We assume that gas (and plasma) has Maxwellian distribution: 




( 8 ) 


and fills the tube homogeneously. We may easily find the kinetic energy and particle fluxes onto 
wall to be: 


q W (n,T)= \ j J v x -^-/ M <fv x dv y dv z = 

— oo — oo 0 

nMVl^Tx 1 ^ 2 

,3/2 J J 


2n p e x p dxdp = 


x ~ o o 

nMVr 3 +°° + r, -r-sj j 

i72~2 ! J S drds = 

^ 0 0 


nkTV T 


n 


1/2 


( 9 ) 


and, respectively. 


+ 00 - 1 - 00+00 

J («.7’)= J J I V X / M <* 

— OO — OO 0 


v x^ v y dv z =\ nV T 


( 10 ) 


Note that this is just kinetic energy flux. Non-mono-atomic neutrals bringing also rotational 
energy to the wall. The expression for total energy deposited is 


w C— 1.5 

Q w = (5 + <r)(a + ) ■ nkTV T 


(ID 


where a=Jt 05 . All particles that collide 
with wall are coming back with wall 
temperature T W =500°C. 


H* 


E, eV 

13.6 


Temporal evolution of the density and 
temperature of the 6 plasma species: 
{ e , H 2 , H, H + , H 2 + , H/} is determined 
by cold gas puff rate, electron heating 
rate, energy exchange with wall, 
mass/energy exhaust and various 
inelastic collisional processes. 
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Figure 17 - Energy levels of H-species 
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Inelastic and elastic collisions . We include the most important inelastic process caused by 
electron impact (electron mobility is highest). The energy structure of the hydrogen species as 
adopted to enforce conservation in present model is presented in Fig. 17. 


Hydrogen molecule dissociation (requires 4.4eV) by electron impact (rate DH2) 

H 2 +<?— ► H+H+e (12) 

<IDH2> = <8.5,11.7,11.7 eV> « lOeV. We assume the electron loses lOeV and each H atom 
carries away H 2 /2 + l/2EDH2(=(10-4.4) = 5.6eV) 

Hydrogen molecule ionization by electron impact (rate IH2) 

H 2 + e t — * e + H 2 + + e (13) 

Electron loses HH2= 15.4 eV, secondary electron is created with OeV energy, H 2 + carries away 
all H 2 energy. 

Molecular hydrogen ion dissociation (takes 2.6eV) by electron impact (rate DH2+) 

H 2 + + e — >e + H + +H (14) 

{From averaged vibrational states 1=2.4 eV, atom and atomic ion carry away extra <E>=4.3eV a 
piece in theory}. We assume electron loses IDH2+=11.2eV, and that atom and ion carry away 
H 2 /2 + l/2EDH2+(=8.6eV) apiece. 

Hydrogen atom ionization by electron impact (rate IH) 

H +e.— > e + H + + e (15) 

Electron loses IIH = 13.6eV, secondary electron at OeV. H + carries H energy. 

Hydrogen atom excitation by electron impact (rate EH) 

H +e. — > e + H*(n=2) — * e + H + y (16) 

Electron loses IEH = 10.2eV, energy ZEH (=10.2 eV) radiated as Lyman-a. 

Heavy ion formation (the only non-electron reaction, rate FH3) 

H 2 + +H 2 — ► H/ + H (17) 

Exothermic reaction with EFH3 = 1.2eV, evenly divided between products. 

We can add the following comment on the included plasma processes. If plasma temperature is 
low at about leV then negative ions H- and vibrational-rotational levels dynamics of H 2 (v,j), 
H2+(v,j) may be important. Other excited levels of H (n=3,4,..) may be important if plasma 
density goes up. Among other elementary processes to be of importance: resonance charge- 
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exchange H 2 -H 2 + and H-H\ 3-body and photo recombination. Coulomb collisions, inter- and 
intra- elastic collisions between neutral particles and ions. Another interesting possibility is H 3 + 
ion formation, but this is non-electron reaction. As plasma density is low ~10 12 cm‘ 3 and 
temperature is moderate Te~6eV, we think that 5 mentioned reactions will be dominating. 

Set of balance equations 

We now can derive equations for plasma and gas density and temperature. Analysis shows that 
instead of temperature it’s better to introduce energy density q=CnT. In this case relative 
equations are simpler and explicitly express energy conservation (while density equation 
preserves mass). By looking on reaction rates [9-10] (Fig. 18) and product energies, we can draw 
the following conclusions: 



1) There will be a large population of Frank-Condon neutrals (via both dissociation reactions). 
These neutrals will be effectively transporting heat to the wall because they are not magnetized 
and because tube radius is much smaller than its length. Exhaust velocity of such hot neutrals 
will be significantly higher than that of the cold neutrals. We’d like to simulate hot component 
more accurately. Therefore, we split atomic neutrals into two populations - cold. He, with 
temperature TdleV and hot, Hh, with T>l-3eV. 

2) Hydrogen atom excitation is more probable than ionization. So, a significant fraction of 
energy can be radiated as Lyman-a. Plasma is optically thin, and all radiated power will be 
absorbed by internal elements (antenna, etc.). Molecular ionization has slightly higher cross- 
section and also energy threshold if compared to atomic hydrogen. This tells that there might be 
a significant fraction of H 2 + ions (and possibly H/) in the discharge, depending on the H/H 2 
ratio. 

The system of equations has to be normalized by choosing an appropriate set of dimensionless 
units (typical length is lm and thermal velocity of leV atom). The following units are adopted: 
[n]=10 13 cm'\ [c]=10' 16 cm 1 2 , [T]=[E]=leV, [L]=lm, [m]=M p (proton mass). We may find other 
units to be [V]=(2eWM) 1/2 , [t]=[L]/[V]=0.7ms, [W]=[T]/[t], [R]=l/[n][t], etc. 

Equation for cold molecule density includes gas puff, molecular dissociation and excitation, 
particle exhaust. Equation for energy density contains related energy sinks/sources and also 
energy conversion at the inner wall. We presented only expression for ionization rate. The rest 
are obtained in a similar way. We keep calculated rates in a tabulated form to eliminate need for 
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their re-calculation during simulation. In these units the final form of the equations for H 2e mass 
and energy density are: 

^H2c M n _ « r » « 

U H2 n e n H2c 1 H2 n e n H2c + 


^H3 n H2+ n H2c yH2c n H2c^T 


A 1 

r “"I 


01 00 * 

^ H2^e ) = J4^v O JJJ2 fM( n e>T e 

n e 0 Inz 

dQfi2c _ Ji_n t° 
ol h1 h2 

-(D H2 n e +Jtf2 n e + ^H3 n H2+ )@H2c 

sA n \/H2c , ifl T W„ j/H2h 
°H2 cQh2c*T + °H2c* n H2h*T 


C H2c (S + (t) 

~4oi ' 2 


SB _ ^H2c<~ ■ 

°H2c ~ T~T » ^H2c - — 


Th2c ~ Qh2c^H2c H H2c 

tH 2c 


V‘ 


88 VW* 


(18) 


The equations for the rest 6 plasma components can be easily obtained using similar approach. In 
the final dimensionless form they read are as follows. For molecular ions: 

2+ _ 


dt - I H2 n e n H2c ~ D H2+ n H2+ n e 
~ F H3 n H2+ n H2 ~ y H2+ n H 2+^ 


Yh2+ ~T 

^ 2+ =I H2 n e Qh 2c (22) 

~( D H2+ n e + F H3 n H2^H2+ 

~ $H 2+ n H 2+^H 2+^S 2+ - ^H2+ n H2+^e 

xA _ C H2+ xB _ 1 
°H 2+ : — . °H 2+~T 


(19) 


T H2+~QH2+ /C H2+ n H2+< C H2+~^ 

v t 2+ =^h 2 + /2 > cP + =ylr~il 
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Equations for cold neutrals has extra particle source at the wall due to hot neutrals interacting 
with the wall: 


dn 


dt 


7T~ I H n p n Hc + YHc n HhYT h 7Hc n Hc^l 


He 

T 


A CT + S 

YHc=-rr- 

4Lo 


yB -1 

Yh,- L 


dQllc T « /~t „ 71 7 /He 

~J t - 1 H n e &Hc d Hc n Hc l Hc V T 

T W U Hh 

+ d Hc n Hh T V T 


3 He = 
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OL 4 L 
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* c 4oL 

Vt H c =^ 
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(20) 


Equations for hot neutral account for doubled rate of hot neutrals production during molecule 
dissociation and their conversion to cold neutrals at the wall. 

dn Hh 

= ~^H n e n Hh + ^H2 n e n H2c + ^H3 n H2 n H2+ 

+ D H2+ n e n H2+ ~ 7 Hh n Hh^T ih 

S + 5(7 
7m - t~z 

4oL 

Jt ' = -I H n e Qm + D H 2 n e (Qh2c +n H2c^DH2^ 
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+ ^-5F H3 (n H2 Q H2 + +h H 2 +Qh 2+ +n H2+ n H2^FH3^ 

~S Hh Q Hh V T HH 


/ Hh 
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dL 4 L Hh 



1 m ~ Qm ICm n Hh’ 



( 21 ) 
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Atomic ion equations account for ionization of hot and cold atomic hydrogen, dissociation of 
molecular ions and ion loss through the nozzle: 


dn 

~dt 


— = I H n e C n Hh + n Hc'> + D H2+ n e n H2+ 


- yH+ n H+ C S + ’■> Yh+~ — 

O ' 5 D H 2+ n e (Qh 2 + + n H 2 + E DH 2 + ) ( 22 ) 

- Jh^hJh+Cs* - Sh*» h J, eg* 

Jf A 1 y'-'f ^ 

d H+-— C H+' d H+-~ 

TH+=QH+ /C H+ n H+’ C H*= 


Heavy ion equations balance H 3 + ion production and their exhaust into the nozzle: 


dn 


hi _ 


dt 


~ ^H3 n H2 n H2+ Yh3 H H3^S 


H 3 


1 

r H 2 - L 


dQ 


~ ^^H3^ H H2 Qh 2+ + 


dt 

h H2+Qh 2 +n H2+ n H2^F2) 

' ^H3 n H3^H3^S^ ~ ^H3 n H3^e 

s?A _ 1 oB _ 1 

0 H3-~j^'~H3' °H3~~ 

v“ 3 =Jt h ,i 3. cF=jrTri 

Th3 = QH3^H3 n H3-> C H3 ~ 2 


( 23 ) 


Finally, the density equation for electron, because of the quasineutrality constraint, reduces to a 
simple form. 

At the same moment energy density equation is the most elaborate. It accounts for all inelastic 
processes included into plasma chemistry kinetic, electron heating and electron energy carried by 
ion sub-sonic flow. 
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I 


« =n H+ +n 


l H2+ + n H 3 + 


d Q e _ w 

dt OL 


' n e ^DH2^H2c n H2c 


~ n e ^IH2^H2c n H2c 

~ n e I DH2 +Dh2+ H H2+ 

~ n e I IH (I Hh n Hh + I He 11 He ) 

“ n e I EH (E Hh n Hh + E Hc n Hc ) 


- ( n H+Cs + + n H2+C S + n H3 Cs > 

S = — • C , c e =- 

e L e 2 

T -Q 1C n 

e e e 


iH2+ 


,H 3 , 


(24) 


For monatomic particles we used following values of specific heat coefficients: C=1.5, for 
diatomic species C=2.5, and for heavy ions C=3. 

Benchmarking against analytical results. For relaxation oscillations see ref. [20]. 


Benchmarking against VX-3 experimental data. 

We performed runs for the following fixed parameters of the helicon deuterium discharge: 
R=5cm, L=70cm, T w =0.05eV, To=0.03eV. We have varied input power about W=l-2.5kW and 
mass of puffed gas, p=50-130sccm. We tried to be close to experimental data: Te=5-7eV, n=l- 
2xl0 12 cm' 3 . Best agreement was achieved under the following assumptions [16]: 

(a) mean gas flow in the source is sub-sonic at M=0.03; 

(b) plasma ‘’parallel” flow velocity VZ=0.3Cs; 

(c) electrons absorb 40% of input power. 



Figure 19 - density vs RF power 



Figure 20 - density vs flow rate 


i 
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In the Fig. 20 plasma density scan versus gas flow rate at fixed input power is presented. Figure 
19 illustrates dependence of the plasma density versus absorbed RF-power. Agreement between 
simulation (solid curve) and experiment is better than 10%. Numerical model explains why 
plasma density goes up faster than the input power. It turns out that electron temperature goes up, 
making the fraction of the energy that goes into gas ionization higher with respect to radiation 
(Ly-a) loses. 

As one can see, the 0-D mass-power balance model reproduces VX-3 experimental data quite 
well. The only difference is observed beyond certain input power per unit mass throughput. In 
the simulation we see that discharge either decays or oscillates. However, in the lab experiment 
this phenomenon was never observed. 


Scaling VASIMR helicon operation to the MW-level power 


When studying low-power VASIMR helicon plasma source we neglected radial plasma 
transport as unimportant. However, at higher power densities we no longer can neglect 
the rising radial loses to the walls. There are two major acceptable models for the cross- 
field transport - classical collisional and anomalous Bohm-like transport: 

1.5 


KJK.~ 


qj 2 t 2 

ce * e 


• 4xl0“ 26 


' », A 


t! 5 b 




= 1.5x10 


-8 


(25) 


where 


n.kT, 1.5 
v ~ A l £ ; ~ 

A -L *•' ,,2 _2 **■// 


■K„ 




and 


r yr 

D b =0.06— -~ £ - = 3.2xl0 4 cm 2 /sec 
eB 

The corresponding heat fluxes are given by expressions: 

kT 

q„,~SKy'(kT,)~SK 1 ^ 


« 0 =5W,D,V,n =Str.D,-i 


respectively. 


(26) 

(27) 

(28) 
(29) 


Estimates show that for 20x1 00cm tube, 4xl03cm-3 plasma density, lOOeV temperature, 
and lOOOG-field strength, the classical flux is negligible, while anomalous 


can’t be ignored.. 


< 7 * =10 


-12 


ST\ 

RB 


= 250.000W 


(30) 


VASIMR helicon source operation simulations indicated possibility of an interesting 
transition to high-performance regime, see Fig. 9. 
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Figure 21 - High-Power scan ofVASIMR discharge parameters 
indicating possible transition to a high-performance regime. 



Figure 22 - Dependence of the VASIMR 
ionization cost versus input power. Also 
indicated is the presently achieved cost. 


Simulations also indicate low, on the order of 60eV ionization cost for such propellant 
burnout regime, as follows from Fig. 22. 
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Conclusions 


Several models have been developed under the reported period to study propellant flow and 
plasma conditions in the helicon plasma source of VASIMR thrusters. 

The models have been benchmarked against VX-10 experiment. Numerical results demonstrate 
good agreement with the lab data. 

The results indicate favorable VASIMR helicon source scaling to MW-level powers. The 
theoretically achievable 


Nomenclature 


B magnetic induction (0-7 Tesla) 

e electron charge ( 1.6 10 19 Coulomb) 

E electric field (Volt / m) 

/ distribution functions 

W heating power ( 0-10 4 W ) 

I sp specific impulse (5000 - Iff* s) 
j particle flux (0 - Iff 21 / m 2 s) 

k Boltzmann constant (1.38 Iff 23 J / K) 
kT/e temperature (1 - 100 electron-Volt) 

L plasma length (1 m) 

c cross-section(m ) 

S quartz tube surface (m ) 

M particle mass (proton mass: 1.67 Iff 27 kg) 
\x propellant flow rate (Iff 6 kg/s) 
n density (0 - 10 19 m' 3 ) 

N density (0 - 10 21 m' 3 ) 

R radius of a boundary (m) 

r L Larmor radius (Iff 2 -7 m) 

t time (s) 

V exhaust velocity (10 4 - 10 5 m/s) 

V T thermal velocity (ion: Iff* - 10 s m/s) 

Cs sound velocity (ion: Iff* - 10 5 m / s) 

aj3Sy coefficients in the equations 

q heat flux (W / m 2 ) 

A d Debye length (W 4 -7 m) 
m fp mean free path (0.7-7 m) 

Q energy density (W / m 2 ) 

Subscripts 
e electron 

/ ion 

H+ hydrogen ion 

H 2 + molecular ion 

H 3 + heavy ion 

H 2c cold molecule 

H c cold atom 

H h hot atom 

p plasma 

RF radio-frequency 
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Conclusions 


Several models have been developed under the reported period to study propellant flow and 
plasma conditions in the helicon plasma source of VASIMR thrusters. 

The models have been benchmarked against VX-10 experiment. Numerical results demonstrate 
good agreement with the lab data. 

The results indicate favorable VASIMR helicon source scaling to MW-level powers. The 
theoretically achievable 


Nomenclature 


B magnetic induction ( 0-1 Tesla) 

e electron charge ( 1.6 10' 19 Coulomb) 

E electric field (Volt / m) 

/ distribution functions 

W heating power (0-10 4 W ) 

l sp specific impulse ( 5000 - l(f s) 

j particle flux (0 - id 11 / m 2 s) 

k Boltzmann constant ( 1.38 10' 23 J / K) 

kT/e temperature (7 - 100 electron-Volt) 

L plasma length (7 m) 

a cross-section(m ) 

S quartz tube surface (m 2 ) 

M particle mass (proton mass: 1.67 10' 27 kg) 
p propellant flow rate (7 O' 6 kg/s) 

n density (0 - 10 19 m' 3 ) 

N density (0- 10 21 m' 3 ) 

R radius of a boundary (m) 

r L Larmor radius (70' 2 -7 m) 

t time (s) 

V exhaust velocity (10 4 - 10 s m/s) 

Vf thermal velocity (ion: 10* - 10 s m / s) 

Cs sound velocity (ion: l(f - 10 5 m / s) 
aflSy coefficients in the equations 
q heat flux (W / m 2 ) 

X D Debye length (10' 4 -7 m) 
mfp mean free path (0.7-7 m) 

Q energy density (W / m 2 ) 

Subscripts 
e electron 

i ion 

H* hydrogen ion 

H 2 * molecular ion 

H 3 * heavy ion 

H 2c cold molecule 

H c cold atom 

H h hot atom 

p plasma 

RF radio-frequency 
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Plasma thrusters are the future 
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MIT 


Nozzle converts T into l„ 
T (eV) * (I /10Y/z 


• Ideally high variable specific impulse 
i sp =E/g is required 


From Saha eq. 

n 2 = 10 22 NT ,s exp(-E/T) 

for H-0 n/N«1% 

! Higher => Increased n/N ■ plasma 






There are -80 types of plasma thrusters MIT 


PIT Isp - 5 Ksec Hall effect thruster Isp - 4 Ksec 

proposed by Morozov in 70s at KIAE 



Ion thruster Isp ~ 3-10+ Ksec 



Magneto-plasma-dynamics -3-10+ Ksec 
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VASIMR: Variable Specific Impulse Magnetoplasma Rocket |Y 



VASIMR was proposed at MIT -1988 

Franklin Chang-Diaz 


electromagnets 


quartz tube 



Helicon source invented in 1968 at Flinders University, South Australia 

Rod Boswell 
















Kinetic study confirms viscous sub-sonic propellant flow MIT 


-V,—+V v — = — — — + Wall + Source 
X dx y dy T(f) 


Eddy formation is possible, non-Maxwellian NDF, but overall M-0.04-0.1 


Kn» 1/4, Re -2 


Kn- 1/8, Re -4 


0-D model of the VASIMR plasma source (H) 


MIT 


Power and mass model of a helicon discharge assumes 

H ; (D ) at room T is puffed at a given mass flow rate Ji. 

a fraction of RF-power P is absorbed by electror s* l * c * OTh *’ 0ns s 

\ f T W" OOSeV 

quasineutrality =>n = n w + nmt + n mt /A A /' 


I \ S 


* plasma is magnetized, doesn’t hit the wall 

- hot neutrals hit wall and reflect with wall T 

- gas and plasma leave the source at M«0.1 & 1 


H 2 infection at 0.025eV 


all 7 species {e, H 2 , H 2 , H c , , H , H 3 } have HjjH 

Maxwellian distributions HkM 

most important plasma-chemistry reactions are included 

> Set of 13(+1) non-linear O.D.E. with coefficients accounting for^^® 
physics and dimensions of the helicon source determines the temporal 
evolution of T and n of each of the plasma species 











Most important plasma chemistry, He 


MIT 


Helium atom ionization 

He + e e + He* + e 
Helium atom excitation 

He + e -» He* He + y 
Ion - neutral atom charge-exchange 
He* + He -» He + He* 
Helium ion ionization 

He* + e -> e + He 2 * + e 
Helium ion excitation 

He* + e He* * -> He* +y 
Wall conversion 

He* > He 



10' t, ev io* 


T has to be > 6eV and <20eV 
' to avoid ion radiation and 
double ion formation 


0-D plasma kinetic model results for D 


MIT 


D discharge for p=50-1 30sccm flow rate and P=1-2.5kW input power has 
shown very good performance. Can our 0-D model reproduce it? 


If (a) we take results from gas flow simulation that gas flow has M=0.06 

(b) assume that plasma ‘’parallel” flow velocity V z =Cs 

(c) guess that 50% of input power is absorbed by electrons 
the agreement with experiment is fare 


© experiment 


500 1000 

P, w 


n vs P (|x fixed) & 


experiment < 



50 60 70 80 90 100 110 120 130 140 150 

D, gasflow, seem 

n p vs \l (P fixed) 









MW-class VASIMR helicon operation, D 


MIT 


Balance model with anomalous diffusion is applied 
to helicon source with following parameters 
0 x L = 0.3 x 1 m, p=0.1g/sec, P-0.1-0.5 MW 


n p x 10 13 cm‘ 3 


fuel eff. 


/ _ energy eff. 

/ Ly/P 

•i l .... i .... i ... I .... ■ 

0 0.1 0.2 0.3 0.4 0.5 0.6 

P, MW 


If RF coupling to plasma is high at high 
powers, then high source performance 
is achievable 


• Plasma density n~ 2 x 10 13 cm 5 

• Propellant efficiency (cold gas to hot 
plasma conversion) > 90% 

• Energy efficiency (gas ionization and 
plasma heating vs heat and radiative 
losses) > 80% 


Single-stage VASIMR is a capable thruster 


MIT 


Operational regime transition occurs when most of the gas is ionized 
|i=0.1g/sec = 1.5 xio” at/sec => 3 KW is required to heat plasma by leV 


T a ,eV 


Thrust, N 


Isp, Ksec 


If RF coupling to plasma is high 
(has vet to be demonstrated at high 
P) then high thrust is achievable 


thrust efficiency . Thrust on the order o{ several N 


0 0.1 0.2 0.3 0.4 0.5 0.6 

P, MW 


• l sp ~ 5-7 Ksec 


• Thrust efficiency (beam energy to 
applied electrical energy) >50% 




/ 

Ionization cost for VASIMR source 

MIT 

n 


Current VASIMR-3kW helicon 
experiments show -300 eV/ion 
production cost for H and He 


eV/ ion 


present VASIMR-3kW 


prediction 


burn-out 


0.2 0.4 0.6 

Power, MW 

Projection: -60 eV/ion for H 

~120eV/ion for He 


Other EP system costs 
Hall Effect Thruster 
SPT-140 '97 -155 eV/ion 
D-100’96-163 
T-160 ‘95-200 
P-5 ’01 - 183 
Ion Thruster 
IAPS ’83 (Hg) -258 eV/ion 
NSTAR ’98 (Xe) — 1 85 
XI PS-25 ’98 (Xe) - 1 15 
Higher for MPD and PPI 


J. of Propulsion and Power , vol. 14, #5, 1998 
and J. Haas Thesis, U. Michigan, 2001 


ASPL data for D: near 100% gas utilization 


MIT 


> The measured ion flux scales linearly with power and gas injection rate below a 
saturation value 

> When the gas flow is below the saturation value the electron temperature rises 
substantially as the flow decreases, indicating a high degree of ionization 

> At the saturation values, the plasma flux equals the input neutral particle flux, 
within error 

> Too much input gas degrades the discharge substantially 



1 ! 




ration 






• saturation •Hm 

i Wmm 


m iSyfUi-v^ < v- : : : 











Why gas burn-out was not seen before? 


MIT 



h 

I 2 

I 3 

Ar 

15.8 

27.6 

40.7 

Kr 

14 

24.4 

37 

Xe 

12.1 

21.2 

32.1 


• Ionization and gas utilization efficiencies are not important for plasma processing uses 


• Most of the experience is with heavy gases 

-Ar, Kr, Xe 

• Radiative loses are dominant, 

full ionization is impossible! 


• Energy and propellant utilization efficiencies are crucial for space applications 

• Hydrogen (D) can be completely ionized. | 

Dark plasma is feasible Li 

• There are other interesting possibilities - M 

u, cs ... h 2 o, nh 3 He 24.6 54.4 


Gas retention by a choke, in conjunction with magnetic field mirror 


H 

13.6 



He 

24.6 

54.4 


Li 

5.4 

76.6 

122 


Gas Choke Installed 


>The same plasma flux with half the 
input gas flow 

>Higher Mach numbers measured 
>Much lower background neutral 
pressures (-10- 5 torr) 
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Predictive VASIMR modeling requires kinetic analysis Pfll I 


VX-10 experiment set-up 

System of 2D2V kinetic equations 


dt [ $y 2 \ m 

Cee + C(i +EX e +l e +H e -R e 



Model geometry 






. v± 3y. 


Cjf + I i + CX i + H i ~R i 

C NN~ I N +CX N + S N 


* 

: i/t\ 


Takes into account Coulomb, various 
elastic & inelsatic collisions, ambipolar and 
sheath potentials, mirror force, heating. 











Effects of ambipolar E-field and B-field trapping 


MIT 


T. Glover, 2001 


=3200s =4200s 




For T =6eV & n m<u 

^ Alj P =1 000 

u min 

s 


Braginskii ambipolar electric field 
along the open magnetic lines: 

e</> = -\.l\VT e -TV\nn p 

<t> « 30-50V — > l jp * 4-6Ksec 


! First stage of VASIMR is an effective thruster 


Conclusions 


MIT 


★ 0-D plasma chemistry model & 1-D hybrid semi-analytlcal mixed collisional gas 
flow models are developed for the VASIMR plasma source 

★ Quantitative analysis of mass/energy balance in VASIMR helicon source is done 
for H,D&He. Good agreement with VX-3 and VX-10 experimental data is achieved 

★ 2D2V purely kinetic model for internal rarified gas flow in the system of channels is 
proposed and verified, shows possibility of eddy formation 

★ Results from kinetic model for plasma and gas show deviation of IDF from 
equilibrium. Its shape is in agreement with the RPA data 

★ Helicon plasma source + extended magnetic field produce ion flow 
with high axial drift velocity due to magnetic mirror and ambipolar potential 

★ Single-stage VASIMR in the VX-1 0 configuration is capable of achieving 
l s/) >4Ksec & thrust -10mN with -50% efficiency. 

★ Model has predicted transition to the gas burnout regime, which has been 
observed later in the ASPL experiments 

★ Model predicts favorable scalability to the MW-level operation with 60eV/ion cost, 
>50% thrust efficiency and >7000s 1^ 








